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The contact value theorem for Coulomb gases in planar or film-like geometries is derived using 
a Hamiltonian field theoretic representation of the system. The case where the film is enclosed 
by a material of different dielectric constant to that of the film is shown to contain an additional 
Casimir-like term which is generated by fluctuations of the electric potential about its mean-field 
value. 



I. INTRODUCTION 

Coulomb gases arise in a huge variety of physical contexts from plasmas to soft condensed matter systems 0, . 
Ideal Coulomb gases, where only point like charges are considered, present a number of useful sum rules 0,3- These 
sum rules are exact identities between certain statistical mechanical observables. Sum rules are useful for checking the 
validity of approximation schemes, which must almost always be applied in the case of strongly interacting systems. 
They also provide potentially useful experimental and numerical methods of indirectly measuring local observables 
in terms of macroscopically measurable quantities such as the pressure and vice-versa. Finally sum rules may also 
be used to verify the accuracy of numerical simulations where, thanks to the sum rule, the same quantity can be 
measured in two independent ways. 

An example of such a sum rule is the contact value theorem which relates the surface charge and density at the 
surface of the system to the system's pressure In this paper we use a field theoretic approach to show how the 

basic contact value theorem can be derived in the case of layered geometries, such as soap films, and show how the 
contact value theorem is modified when the dielectric constant (e) within the film is different to that outside the film 
(eo). The film geometry is of particular importance as it corresponds to the experimental set up used to study soap 
films [j| . We exploit the planar geometry to develop a Hamiltonian formulation of the Sine-Gordon field theory which 
arises for Coulomb gases In this formulation the perpendicular direction, denoted here by z, acts as a temporal 
coordinate in which a field 4>{r) on the plane perpendicular to z, and whose coordinates are denoted by the vector 
r, propagates. The case where global electro-neutrality holds is treated within this formulation. The condition of 
electro-neutrality can be related to the choice of ground-state wave functional for the dynamical field </>(r). 

In the canonical ensemble, the partition function for a system of fixed particle number on a 3-d space V is given by 
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Tr J d[(f)}exp(-~ J dx e(x) f Vc/>(x) \ 2 + i/3 J dx p 



Z = Tr / d[4>] exp ( -'- / dx e(x) ( V0(x) ) 2 + i(3 / dx p e (x)^»(x) ) , ( 1 ) 

where Tr denotes the classical trace over the particle positions -^j J v YiiLi an d d[4>] denotes the functional integral 
over the field (j>. The field <f> is the Wick rotated electrostatic field yb and is thus related to ip via yb = —i<f>- The 
term p e ( x ) = J2iLi ~ x i) + P<j( x ) is the charge density of the system. The first term is the dynamical charge 
density which can vary in the system, xj being the position of particle i and qi its charge. The second term p q 
is a quenched background charge which is not dynamical and represents, for example, a fixed surface charge. The 
integration volume in the action V is all space. We note that the above treatment of a two component Coulomb 
gas needs to be modified where the basic physical description of point charges interacting via a Coulomb potential 
is thermodynamically unstable. For instance the system can become unstable and have a tendency to collapse at 
low temperatures if some short range, for example hard core, repulsion is not included. The Sine-Gordon theory can 
always be regularized by introducing a high momentum or short distance cut-off in the Fourier modes of the field 4>. 
We note that the above formulation contains self interactions between the particles i.e. the terms qiqjv(r{ — rj)/2 for 
i = j where v is the effective pairwise interaction are included. The interaction of a particle with its image charges 
is part of this contribution and should be included, but the self-energy in the bulk medium should not contribute to 
the physical pressure and so must be subtracted. For a monovalent system in dimension d if one removes the bulk 
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self-energy term the Sine-Gordon free energy is corrected by a term 



AF = -iY^ = -^A/W 3 (2) 



2e (2nYJ Q 

where N is the number of particles, Sd the surface of the unit sphere in d dimensions and A is an ultra-violet or short 
distance cut-off. Note that for d > 2 AF is a regular function of A and so can be absorbed into the fugacity. However, 
for d — 2 the integral in Eq. J5J) has an infra-red divergence at k — which must be cut-off at the inverse system size 
l/L giving 

AF = -^(ln(A)-ln(l/i)) . (3) 
Thus, in two dimensions there is a correction to the Sine-Gordon pressure of 

AP = -^ (4) 

where 2p = N/L 2 is twice the density of electrolyte.. 

Normally the particles i are restricted to a sub-volume V, for instance in the interior of a soap film for electrolyte 
solutions. In the case where e is constant, the functional integral over </> is easily done and we recover a system of N 
particles interacting via the Coulomb potential. For a varying dielectric constant e the resulting interaction depends 
on the spatial variation of e and the resulting integration gives rises to a more complicated interaction which can 
be interpreted in terms of image charges. The partition function Z in Eq. thus contains a term which is due 
to the pairwise inter-particle interaction, plus a functional determinant coming from the integration over <j). Both of 
these terms are present in the physics of the problem and should be taken into account. The form of Eq. (JJJ comes 
directly from a static approximation to Quantum Electrodynamics (QED) where the non-zero frequency Matsubara 
frequencies and electric currents (and thus the magnetic field) are neglected in the action of QED |E S 13 ■ This 
approximation is justified when the charge distribution is very weakly be coupled to non-zero Matsubara frequencies 
0, that is to say the response time of the charge to the non-zero frequencies is large. The non-zero Matsubara 
frequencies, when decoupled from the charge distribution yield a van der Waals interaction between the surfaces in 
the problem and can be calculated independently in this approximation . 



II. HAMILTONIAN FIELD THEORETIC FORMALISM 



For simplicity in what follows we consider a film system, that is to say an inner region having two interfaces 
separating it from an outer region. The surface of the film has has area A in the (x, y) — r plane, the direction z is 
perpendicular to the film surface. For a system consisting of a film of thickness L, the film is in the region z € [0, L] 
and the exterior of the film is in [— T',0] (the left exterior) and [L,T — L] (right exterior). The total length of the 
physical system in the z direction T + T" is taken to be constant and we consider the limit of T and T' large. In 
the simplest case, which we study here, the electrolyte is monovalent and the fugacity of cations and anions is the 
same and denoted by \i within the film region and is zero outside. This can be encoded in a spatially varying fugacity 
/x(x) = /j,(r,z) = fJ-(z), with /i(z) = \x for z € [0,L] and fi(z) — for z (fc [0,L]. The dielectric constant within the 
film is denoted by e and the external dielectric constant is denoted by e . For a soap film, for example, eo could be 
the dielectric constant of air and e the dielectric constant of water. We may also consider systems where the regions 
z G [—(5, 0] and z S [L, L + 5] have a dielectric constant determined by the dielectric constant of a surfactant and its 
concentration at the surface. 

In the grand canonical ensemble the grand partition function is written as 



The film the pressure is therefore given by 



d[0]exp(-5[0,L]). (5) 



In the case considered here, S is the action of a generalized Sine-Gordon field theory 6 

S[<p,L] = ^ [ dx e(z)(V(j}) 2 - 2 f dx /i(z)cos (e/30(x)) 
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~ i/3 I dr a<l>(x) + ifi I dr o-0(x), (7) 

Jz=0 Jz=L 
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where the integrations over 3-space, denoted by dx = drdz, are over the (x,y) = r plane of area A and over the 
coordinate z. The last terms are the (constant uniform) surface charge contributions from the transverse planes of 
area A at z = and z = L. 

We now rewrite the action using the 2 + 1 decomposition in terms of the field (f>(r) which evolves with a temporal 
coordinate z 



E = J V[4>] exp (- J dz S[(j), z] 



(8) 



where the path integral action S is given by 



dr- 



pe(z) 
2 J A \dz 

i(3a (5(z)+S(z - L)) 
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dr o 



[ dr (V r 

J A 



2/i(z) / dr cos (e 

I A 



The functional Schrodinger Hamiltonian for the path integral outside the film is given by 

1 S 2 



Inside the film the Hamiltonian is 

H F = I dr 



dr 
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2(3e 5(j}{rf 



4>) — 2/icos (e[34>) 



(9) 



(10) 



(11) 



The terms containing the surface charge and other more general surface interactions may be included as source terms. 
We note here that in the case of the 1-d Coulomb gas the field 4> is interpreted as the position of a single particle and 
the quantum mechanical formalism then leads to an exact solution [10| . In this formulation we postulate that 



5 = Tr exp {-TH E ) O exp {-LH F )) O exp (— (T' - L)H E ) , 



(12) 



where here Tr denotes the trace over a complete set of states. The source term O for a constant surface charge is 
clearly given by 



O = exp ^if3a J dr 0(r) 



(13) 



The reason we say that Eq. I|12l) is postulated, is that the boundary conditions for the above path integral are not 
straightforward to determine. We shall see later that the above choice assures the electro-neutrality of the film since 
in the limit T, T' — > oo it only involves the ground state wave functional |^o) of t ne Hamiltonian He- In general, 
instead of taking the trace, we could specify any linear combination of wave functionals of He as the initial (and final) 
state for the field <fi, however as long as it has a non-zero component of the wave function |Wo), in the limit where T 
and T' are large the result will be the same. So, in the limit of large T and T' the grand partition function may thus 
be written as 



5 = (* | exp (-TH E ) Oexp (—LHf) O exp (— (T" - L)H E ) |*o>» 
In this formulation, the derivative with respect to L is easily and unambiguously taken, we find that 

|= = -(* | exp {-TH E )0 exp {-LH F )) H F O exp {-(T' - L)H E )\V ) 
+ (* | exp {-TH E ) O exp (-LH F ) OH E exp {—{T' — L)H E ) |*o). 
We now define the momentum operator for the field (f> at the point r by 

. S 



(14) 



Pa, 



(15) 



(16) 



This leads to the commutation relation 



[P <Kr) ,<j>(r')}=-i8(r-r'). 



(17) 
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The kinetic operator K is then defined by 



K =A* p *«=-//-w (18) 



In this notation the exterior and interior Hamiltonians read 



2/3e 

where the functional potentials of the Hamiltonians are given by 

(3e 



HE = Wo +VE (19) 
Hp = ^- + V F , (20) 



V * = 2 

v - Pe 
V F - - 



[ dr (V r 0) 2 (21) 

J A 

I dr (V r </>) 2 - 2/i / dr cos (e(3<f)(r)) . (22) 

J A J A 



We note that these functional potentials are pure functionals and involve no functional derivative operators, thus they 
commute with other functionals, notably O. The Eq. Ijl5|l is now written as 



dL 

Using the result 



(*o| exp (-TH E ) Oexp {—LHp) ([Hp, O] - 0{H E - Hp)) exp (— (T' - L)H E ) |* >. (23) 
[P 2 } , O] = 2(3aOP 4>(r) + (3 2 a 2 0, (24) 



we obtain 



8Z A (ia 2 
dL 



-A^-E - ^<* | exp (-TH E ) O exp {-LHp) O Qf drP^ exp (— (T' - L)H E ) |¥ ) 
+ (* | exp (-TH E ) O exp {-LHp) 0(H E -H F ) exp {-{T 1 - L)H E ) | *„). (25) 
The final result for the pressure is thus 

= + ^<(#E - H F )\ S+ ) -J^iJ/r f«r)U), (26) 

where the above notation indicates the normalized expectation values of the operators shown, evaluated at the 
rightmost outer-surface of the film s + i.e. at z — L + . The third term above can be shown to be zero in the case of 
systems which are globally electro-neutral; we will demonstrate this more technical point later. We note the relation 



Hp = —H E - dr 
e J A 



(3e / e 



2/icos (e/30) - y ( 1 



I) (Vr 



(27) 



Also it is straightforward to see that 

( M exp(±ze/?0(r)| z ) Hp±(r,z)), (28) 
is the mean value of the cation/anion density at the point (r, z). Putting these results together we obtain 

(3P = {p+\ s + p -\,)-P?l-t(i-&) ((V r 0) 2 | s ) + i (l - C j) (H E \ S+ ), (29) 

where all but the last term are evaluated on the surface s (z = L) at any given point (by homogeneity in the plane 
A). Here there is no ambiguity with regards to the interior or exterior of the surface as the terms are pure functionals 
of the field <f> and commute with O as it is. itself, a pure functional of the field <f>. The last term is evaluated at the 
outer-surface s + (at z = L + ). We thus obtain 



(30) 
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where Eq is the energy of the ground-state wave functional |\&o)- ln the case where eo = e, Eq. I|30|) immediately 
yields the classic contact value theorem as the third and fourth terms are identically zero. 

While the functional terms appearing in the expectation values of the Hamiltonians are easy to interpret in terms 
of observables, the kinetic term K requires more work. The key result here is 



' 1 p2 I 

v 2/?e(z) *«' 



x 2 \ dz 



This can be seen in the Heisenberg formalism which gives 



dz 

We therefore have that for a point z inside the film 

/3e / d(j)(r, z) 



(HfU 



dr 



2 V dz 



(f>) — 2^cos (e(3(j>) 



and outside the film 



(H 



E z = 



dr 



/3fo ( d<f>(r, z) 
2 



dz 



Pe 



(V r 



Using Eq. I|34|) in Eq. Ij29(l we find the alternative expression 



/3P= ( P + \ s +p-\s) 



/3a 2 
~2~T ~ 2 



( e -e )((V r ^) 2 | s ) 



(3e 



| s+ >- 



(31) 



(32) 



(33) 



(34) 



(35) 



In mean field theory, when the field <fi is replaced by its mean field electrostatic field —i4>c, the electro-neutrality of 
the system implies that dip c /dz = outside the film and the homogeneous nature of the mean field solution in the 
plane of the film yields V r -0 C = 0. Thus the third and fourth terms of Eq. (|35l) are zero and hence the contact value 
theorem as classically stated is always verified at the mean-field level. If we expand about the mean field solution, we 
see that the correction term due to the variation of the dielectric constants comes from fluctuations about the mean- 
field solution. The value of this term can be calculated in the case of weak electrolyte strength in the Debye-Hiickel 
approximation [^li^ll"l||. 

We now return to the question of the global electro-neutrality of the system. Using the correspondence between 
the field <fi and ip we have 



dz 



dz ' 



(36) 



If the film is electro-neutral then, by symmetry about the film's mid-plane, the integral of the electric field over every 
plane perpendicular to the z direction outside the film must vanish, and thus 



dr 



dipjr) , 
dz 1 



= 0, 



for z > L (and also for z < 0). Using the Heisenberg formalism we find 

dcj)(r) 



dr 

A dz 



zij^dr [H E ,<j>(v)}\ z ) = -L(^drP(r)| z ), 



and using the fact that 



cxp (— (T' - L)H E ) |*o> = exp (— (T' - L)E ) |* ), 



we may write 



f drP m exp(-(T'-L)H E )\y )=eiq>(-(T'-L)H E ) f drP 0(r) |* o ). 

J A J A 



(37) 



(38) 



(39) 



(40) 
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The final step in the proof of our result is to show that 

dr P 0(r) |* o > =0. (41) 

Outside the film we consider the Fourier representation of the field 

= 4=^<Kp)cxp(* P t). (42) 
* p 

In this basis the exterior Hamiltonian becomes 

2 

p 



p)<M-p) 



(43) 



Thus outside the film each Fourier mode of <p is the coordinate of an independent Harmonic oscillator and the ground 
state wave functional is given by 



(0|* o ) <* II ex P (-^ MW (P)^(P)^(-P)J . ( 44 ) 
P \ / 

where M = (3eo and w(p) = |p|. We now note that 

drP 0(r) = -n/I-^—. (45) 

00(0) 

However from Eq. I|44(l |^o) is clearly independent of 0(0) and so the desired result Eq. (|41|) follows directly. 
Using the fact that the Hamiltonian H F commutes with itself we may also rewrite Eq. ()15f) as 

|= = -(^ \cxp(-TH E )Oexp ^LH F )^ H F exp ^LH F )^j Oexp(~(T' ^ L)H E )\^> ) 

+ (^ \exp(~TH E )Oexp(-LH F )OH E exp(-(T' -L)H E )\*o). (46) 
Applying the results derived above and using Eq. leads to 

pp = ( P+ \ m + p-U f <(v r 0) 2 | m ) + f <(g ) 2 L> + ^((v r 0) 2 1.) - ^<(f£f U>, (47) 

where to in the above indicates the value taken at the mid-plane of the film z = L/2. At the mean field level for 
an electro-neutral film the last term again vanishes and, in addition, by symmetry the mean field solution tp c obeys 
dijjc/dz — at z = L/2. Hence if we neglect the corrections to mean field in the second two terms, the so called 
mid-plane formula for the mean-field approximation 0, Q is immediately recovered. It is also worth noting that Eq. 
(|47|l is valid for any surface charge operator O, thus explaining why various charge regulated models obey the classic 
mid-plane formula at the mean- field level [y, [lj • 

The Hamiltonian formalism is also illuminating in the case of a bulk calculation. Here the grand partition function 
is given by 

E B = Trexp(-LH F ) (48) 

which gives the pressure bulk pressure as 

0Pb = -j{H f ) (49) 

where the expectation above is taken at any value of z. Now using the relation Eq. ll.''>.'iH and the homogeneity in the 
transverse plane A we find 

p PB = 2 P -^A?) + P idfX) (50) 
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where p is the bulk electrolyte density. Using the isotropy of the bulk system we finally obtain 

/3P B = 2p-^((V0) 2 ), (51) 
6 

where the gradient above is the full three dimensional gradient. The first term above can be interpreted as an osmotic 
term and the second term has its origin in the Maxwell stress tensor for the electrostatic field [13| . 

In d dimensions, using exactly the same decomposition in terms of a temporal coordinate z and ad-l dimensional 
hyper-surface, the above expression becomes 

P p B = Zp-Yd {d ~ mv<t>)2) - (52) 

Note that the above expression recovers 14] the physical pressure for a two-dimensional neutral Coulomb gas which 
is given by the simple formula 

P p = 2pk B T(l - ^-^) (53) 

where the Sine-Gordon pressure has been corrected by the 2 dimensional self energy term of Eq. This result can 
be simply understood from the logarithmic nature of the Coulomb potential in two dimensions and is valid in the 
region of the thermodynamic stability of a purely Coulomb system in two dimensions |l4j |. 

The expression Eq. (|51|l may also be used to calculate the physical bulk pressure of a 3-dimensional Coulomb gas 
within the Debye-Hiickel approximation. One finds that 



where A is again an ultra-violet cut-off and m = y/2pe 2 (3/e is the Debye mass. This gives for large A 

1_ /A3 
12tt 2 It ' 2 



pp B =<2p- — (— -m 2 A + m 3 -) (55) 



which is clearly divergent. However the physical pressure due to the ions P p is given by 

/3P P = PP B - PPb (0) - (3P S (56) 

where (3Pb (0) is the pressure due to an electrostatic field in the absence of ions and Ps is a contribution coming from 
the self interaction of the ions with themselves. To see the origin of Ps we note that 

eV0 = -ieVip = iE (57) 

where E is the electric field. Clearly E = J2i where E^ is the field due to the single particle i. The term in E 2 
which contributes to the interaction between particles is thus E 2 — J^. E 2 . This means that 



i 

Hence if is the electrostatic potential generated by a single particle at the position of interest then 



Pe2f3 f dk = -^m 2 K (59) 



6?r 2 e J " 12tt 
Putting all this together gives us the well known Debye pressure formula 

TTi 

We remark here that the general result Eq. I|52|l may also be derived by putting the field theory on a lattice and 
changing the volume of the system by varying the lattice size |l5j . 



III. CONCLUSION 



We have shown in the case of planar geometries that the Sine-Gordon type field theory for Coulomb systems can be 
expressed in terms of a 2 + 1 path integral for a field cf> on the 2-d surface parallel to the film surface. The expression 
has the form of an S-matrix and the choice of the ground-state wave functional for the ingoing and outgoing states of 
the field 4> imposes the global electro-neutrality of the system. In the case of constant dielectric constant throughout 
the system, the classic contact value theorem is recovered. In the presence of dielectric variations, we see how the 
contact value theorem is corrected by fluctuations about the mean-field solution. A corresponding version of the 
mid-plane formula is also derived, which is also shown to be generically valid at the mean-field level. We again see 
that the corrections come from fluctuations about the mean-field solution. In both cases these fluctuations about 
the mean-field solution are known to be of Casimir type and can be calculated within various approximations. The 
results show clearly that the pressure of the system has two distinct, though inter-related, contributions: an osmotic 
pressure term plus a term coming from static thermal fluctuations of the field 4>. The present Hamiltonian technique 
has the advantage of giving an unambiguous way of taking the derivative of the grand partition function with respect 
to the film thickness L in the presence of surface charges and dielectric discontinuities. In the usual field theoretic 
formulation the taking of this derivative and the interpretation of the resulting terms as thermodynamic averages is far 
from obvious. In this Hamiltonian approach, the passage between the contact value result and mid-plane result is also 
straightforward. We have also shown how the Hamiltonian approach provides an alternative method for representing 
the pressure of bulk systems, giving results in accordance with those obtained via other methods, 

As far as future work is concerned, using the results obtained here similar results can be obtained for more com- 
plicated layer geometries (additional dielectric layers for instance) and also for more complicated surface charges, for 
instance modulated surface charges and surface charges built up from thermodynamic or chemical surface charging 
mechanisms (i.e. charge regulated models) [ftf!^. 



[1] J. Israelachvili, Intermolecular and Surface forces, (Academic Press) (1992) 

[2] W.B. Russel, D.A. Saville and W.R. Schowalter, Colloidal Dispersions (Cambridge University Press, Cambridge) 
[3] D. Henderson and L. Blum, J. Chem. Phy s 69, 5441 (1978); H. Wennerstrom, B. Jonsson and P. Linse, J. Chem. Phys 
76, 4665 (1982). 

[4] M. Desern and C. Holm in Electrostatic Effects in Soft Matter and Biophysics eds. C. Holm, P. Kekicheff, and R. Podgornik, 

NATO Science Series II - Mathematics, Physics and Chemistry, Vol. 46, Kluwer, Dordrecht (2001) 
[5] K. Mysels and M.N. Jones, Discuss. Faraday Soc. 42, 42 (1966); D. Exerowa, D. Kolarov and Khr. Khristov, Colloids 

Surfs. 22, 171 (1987); D.S. Dean and D. Sentenac, Europhys. Lett. 38, 9, 645, (1997); D. Sentenac and D.S. Dean, J. 

Colloid Interface Sci. 196, 35 (1997); V. Castelletto et ai, Phys. Rev. Lett. 90, 048302 (2003J 

[6] D.S. Dean and R.R. Horgan, Phys. Rev. E 65, 061603 (2002);D.S. Dean and R.R. Horgan, eond ma! .0210 195 
[7] R. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill) (1965) 
[8] D. Bailin and A. Love, Introduction to Gauge Field Theory (Adam Hilger, Bristol and Boston) (1986) 
[9] J. Mahanty and B. Ninham, Dispersion Forces (Academic, London), (1976) 
[10] S. Edwards and A. Lenard, J. Math. Phys. 3, 778, (1962); D.S. Dean, R.R. Horgan and D. Sentenac, J. Stat. Phys. 90, 

899 (1998) 

[11] P. Attard, D.J. Mitchell and B.W. Ninham, J. Chem. Phys. 88, 4987 (1988); P. Attard, D.J. Mitchell and B.W. Ninham, 

J. Chem. Phys. 89, 4358 (1988); P. Attard, R. Kjellander and D.J. Mitchell and B. Jonsson J. Chem. Phys. 89, 1664 

(1988); R.R. Netz, Eur. Phys. J. E 3, 131 (2000) 
[12] B.W. Ninham and A. Parsegian , J. Theor. Biol 31, 405, 1971, ;0. Spalla and L. Belloni, J. Chem. Phys. 95, 7689 (1991); 

O. Spalla and L. Belloni, Phys. Rev. Lett. 74, 2515 (1995);M.N. Tamashoro and P. Pincus, Phys. Rev. E 60, 6549 (1999); 

P.A. Pincus and S.A. Safran, Europhys. Lett. 42, 103 (1998) 
[13] L.D. Landau and E.M. Lifshitz, Electrodynamics of Continuous Media Vol 8 Course of Theoretical Physics (Pergamon 

Press) (1960); B. Jancovici, J. Stat. Phys. 99, 1281 (2000) 
[14] A. Salzberg and S. Prager, J. Chem. Phys. 38, 2587 (1962); E.H. Hauge and PC. Hemmer, Phys. Norvegica 5, 209 (1971); 

L. Samaj and I. Travenec, J. Stat. Phys., 101 (2000) 
[15] D.S. Dean and R.R. Horgan, in preparation 



8 



